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Abstract 

In recent years, the Gribov-Zwanziger action was refined by taking into account certain dimension 
2 condensates. In this fashion, one succeeded in bringing the gluon and the ghost propagator obtained 
from the GZ model in qualitative and quantitative agreement with the lattice data. In this paper, we shall 
elaborate further on this aspect. First, we shall show that more dimension 2 condensates can be taken 
into account than considered so far and, in addition, we shall give firm evidence that these condensates 
are in fact present by discussing the effective potential. It follows thus that the Gribov-Zwanziger action 
dynamically transforms itself into the refined version, thereby showing that the continuum nonperturba- 
tive Landau gauge fixing, as implemented by the Gribov-Zwanziger approach, is consistent with lattice 
simulations. 

1 Introduction 

The infrared behavior of the gluon and ghost propagator has received a lot of interest in recent years, 
in particular in the Landau gauge. Many of the discussions were evolved around the zero momentum 
value of the gluon propagator and the infrared enhancement of the ghost. The common belief is now 
that in 4D and 3D the ghost propagator displays no enhanced behavior, while the gluon propagator ex- 
hibits positivity violation, being suppressed in the infrared. Moreover, it attains a non-vanishing value 
at zero momentum. These results are supported by many lattice data [1, 2, 3, 4, 5, 6, 7, 8] as by many 
analytical approaches [9, 10, 11, 12, 13, 14, 15, 16, 17]. Such propagators have been used to extract 
results on the spectrum of gauge theories, see e.g. [18, 19]. In particular, in the Gribov-Zwanziger (GZ) 
framework, which accounts for the existence of (most of) the Gribov copies in the path integral [20, 21], 
this behavior of the ghost and gluon propagator was explained by taking into account the existence of a 
certain BRST invariant dimension 2 condensate [22, 23]. This was called the refined Gribov-Zwanziger 
framework. This particular condensate was investigated as it corresponds to a BRST invariant operator. 
However, one could go one step further. The Gribov-Zwanziger action has a softly broken BRST symme- 
try [20, 22]. Despite this, it is still renormalizable thanks to a wide set of Ward identities obeyed by the 
GZ action. Therefore, one could ask why one would only investigate d = 2 BRST invariant condensates? 

In fact, there exists a whole range of d = 2 condensates overlooked so far, which might be taken into 
account. In this paper, we shall firstly explore these condensates and show that they affect the gluon 
and the ghost propagator, although not altering their qualitative behavior. The gluon propagator is still 
suppressed and non-zero at zero momentum, and the ghost propagator is not enhanced. Secondly, , we 
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shall also be able, for the first time, to calculate the effective action with the help of the local composite 
operator (LCO) formalism at lowest order and give arguments that there is in fact condensation. We shall 
show that the minimum of the effective potential including the condensates is a non trivial minimtmi, 
i.e. in this minimum the condensates are present, leading to a dynamical transformation of the GZ action 
into the refined GZ action. 

This paper is organized as follows. In section 2, we shall briefly review the construction of the Gribov- 
Zwanziger action. The first main point of this paper shall be proven in section 3, i.e. there can be more 
d = 2 condensates affecting the GZ action than considered so far. The second main point of this paper is 
presented in section 4, namely: the construction of the effective action with the help of the local compos- 
ite opeator (LCO) formalism [24, 25]. We first explain the LCO formalism and then apply it to the GZ 
action with the inclusion of the set of d = 2 condensates. We then show that searching for extrema of the 
effective action automatically leads to nonvanishing condensates, i.e. to the refining of the GZ action. In 
section 5, we present the form of the gluon and the ghost propagator and show that they are in qualitative 
agreement with the current lattice data, irrespective of the details of the condensation. In section 6 we 
collect our conclusion. Technical details are provided in a series of appendices. 



2 Summary of the Gribov-Zwanziger formalism 

The Gribov-Zwanziger action takes into account the existence of Gribov copies by restricting the domain 
of integration in the functional integral to the Gribov region Ci, which is defined as the set of field 
configurations fulfilling the Landau gauge condition and for which the Faddeev-Popov operator, 

5W''* = -9p(3;,5"^+g/.icA^) , (1) 

is strictly positive. In [30] it has been firstly shown that this restriction to the Gribov region H can be 
established by considering the following (local) action 

5gz = So-{-Sy (2) 

with 

with 5ym the classical Yang-Mills action and 5gf the Landau gauge fixing 
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5gf = I d^x (b^d^^^ + ^d^Dfc") . (4) 

The fields (^fl^ are a pair of complex conjugate bosonic fields, while (o5^^,(b^') are anticommuting 

fields. We recall that we can simplify the notation of the additional fields ^(p^^, (p][5'^,a5°^, co°'^j in 5o as 
5o displays a symmetry with respect to the composite index ; = {fi,c). Therefore, we can set 

(f-,9-,ra-,a)-) = ((p?,(pf,®f,a)f) , (5) 



and thus 



So = Sym+ % + / d^'x (yld^ [of^^f) - ^d, [of^) - gf^'d^mt , (6) 
The BRST variations of all the fields are given by. 



M« = -(Z)^c)", sc" = -gf''^c''c', 



1 

r- 

^c^=^>^ sb'' = 0, 

smf=lpf, s^f = 0. (7) 
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The massive parameter y, called the Gribov parameter, is not an independent parameter of the theory, 
being determined in a self-consistent way by the following gap equation, commonly known as the horizon 
condition. 



{gf^'A^^^'P + (gr'^AX'} +2y'd{N^ - 1) = 
which ensures the restriction to the Gribov region. This gap equation can also be written as 



ar 



0, 



with r the quantum action defined as 



J [d4>]e 



(8) 



(9) 



(10) 



where /[d$] stands for the integration over all the fields. The action 5gz is renormalizable. For the 

benefit of the reader, we have presented the full algebraic proof of the renormalization of this action in 
the Appendix A, since we have to built on this anyway later on. Let us also mention that, recently, an 
alternative approach was worked out to study the renormalizability of the GZ action [26, 27]. In this 
paper, we shall however follow the original approach of e.g. [28]. 

We recall that the GZ action breaks the BRST symmetry explicitly [20, 22]. This is due to the y-dependent 
term, Sy, and one can easily check from (7) and (2) that, 

^5gz = s{So + Sy) = s{Sy) = -dj d (a>^^ - {Dirn [il' + (1 1) 

3 Further refining of the Gribov-Zwanziger action 
3.1 Introduction 

So far, the GZ action has been refined [22] by investigating the BRST invariant d = 2 condensate 
(^fcpf — mfcoj') and the well known condensate (Apk'j^). The first condensate assures that the gluon 
propagator is non-zero at zero momentum [22], while the second condensate is indispensable in order 
to find a good quantitative agreement with the lattice data, see [4, 29]. The resulting action, called the 
Refined Gribov-Zwanziger action (RGZ), gives rise to a ghost propagator which behaves like l/p^ for 
small p^, and to the tree level gluon propagator given by 



„2 I ■„2 I 2g^N-f 



PliPv 



gafe 



PmPv 



gab 



(12) 



whereby is the mass related to the condensate ((pf (pf — of (Of ) and rrp- to (A^^). We clearly observe 
that this propagator is non-vanishing at zero momentum due to the presence of the mass M^. 



However, as the GZ action breaks the BRST symmetry anyhow, see expression (11), there is a priori 
no need to keep the operators (pf (p" and 55° CO? in a BRST invariant combination, i.e. (^f (pf — (3faff) = 
s (rof cpf ) . In fact, we can spUt the operator into two separate operators, coupled to different sources. 

Moreover, there are also other d — 2 operators, which were overlooked so far. In fact, all possible renor- 
malizable d = 2 operators O; in the GZ action, which have ghost number zero, are given by' 



(13) 



We shall only investigate condensates which are fully contracted over the indices (a,;), e.g. like 9f'(pf = 
'^fi^'ii- However, it is possible to make different contractions over the color indices as is shown in [31]. 

'We are not considering the operator ft" here. A (c"c") condensate would result in massive ghosts, something which is clearly 
excluded by lattice simulations. If iJ^c" is not directly coupled to the theory, it can neither radiatively appear due to a shift symmetry 
of the underlying action, viz. c" — >^ c° + cte, with cte a constant Grassmann parameter. 
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Therefore, if one wants to be absolutely complete, one would have to take into account all possible color 
contractions. Unfortunately, this would be hopelessly complicated. Though, we hope that a good de- 
scription of the IR behavior of the gluon and ghost propagator has been captured by taking into account 
only one color combination. Comparison with lattice data in 3D and 4D seems to confirm this, at least 
so far, [4, 29]. 

We also wish to point out that by including the possibility of condensation of certain operators, we are 
looking at the GZ dynamics w.r.t. a dynamically improved vacuum, in particular an improved calculation 
of the effective action, and thus of the horizon condition via (9), becomes possible. 

3.2 The action with inclusion of J = 2 condensates 

We propose to study the following extended action, 

^CGZ = 5*07 + -^A^ + 'S<p(p + 5j5(o + Sjpp^joip + S<p(p_(o(p + 5vac (14) 

whereby Sqz is given by equation (2) and 



5^2 = Id-^xQA^^-fx- 



L2 
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S<m = j d^xs{V^(ifi) = y"d4;c[lVro>°-V(pfcof] , 

-J d'^x[a{QQ + QW) + P{QW + WW) + xQi + ?,Wx] . (15) 
We have introduced a source t and 4 new doublets of sources, i.e. 





= , 




sP 


= Q, 


sQ = 0, 


sV 


= w, 


sW = 


sG'j 


= iW , 


sW - 


sH'j 


= G'j, 


sG'j = 



(16) 

whereby T is a bosonic source and P, V, H'^ and H'^ are Grassmann quantities. For consistency, the 
sources with double index are symmetric in these indices. In this light, we use the following definition 
for the derivative w.r.t. a symmetric source A^: 

^ = ^(5,.5,, + 5,5,,). (17) 

Notice that some sources have double indices, e.g. H'J , while other sources have no indices, e.g. P. The 
reason for this is only related to the algebraic proof of the renormalization in order to keep certain sjm- 
metries, and has no further meaning. 

We have also introduced a vacuum term, S^^c, which shall be important for the renormalization of the 
vacuum energy. As shown in [24, 25], the dimensionless LCO parameters a, P, %, 5 and C, of the quadratic 
terms in the sources are needed to account for the divergences present in the correlation functions like 

{Oi{k)Oj{—k)), with Oi one of the operators given in expression (13). 

Now we can prove that the action (14) is renormalizable to all orders. The proof is very similar to 
that of the renormalizability of the GZ action, the only difficulty is that the mixing between different 
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sources and parameters is now allowed. We refer to the appendices B and C for all the details. 



For the rest of the work, we are only interested in a restricted number of condensates. Therefore, we 
first set the source W = 0, which is coupled to MCO, as this is not of our current interest^, and we also set 
P = V = r| = 0, as we have introduced these sources only to preserve the BRST symmetry. Secondly, we 
also take H'j = i7'^ = and we set G'^ = 6''G and G'' = WG. The action (14) becomes. 



2cGZ 



5gz + 



7 



1 

" 2 
pGG 



(18) 



whereby {Kd{N^ - \ ) + 'kcfl{N^ - \ f ) was replaced by one parameter p. 



3.3 A diagrammatical look at the potential mixing and at the vacuum 
divergences 

Before starting the calculation of the effective action, we can provide some simplification with the help 
of a diagrammatical argument. Firstly, looking at the action (18), we see that a term xQ'Z is present. 
This term is responsible for killing the divergences in the vacuum correlators (A^(x)(p(p(y)) for x ^ y. 
However, we can prove that there are no divergences of this kind in the one loop diagrams. Let us start 
by considering these one loop diagrams. There is only one possible type of diagram for (A^(x)(p(p(>')), 
as displayed in Figure 1. 




Figure 1: 1-loop diagram for (A^(j[:)(p(p(y)). 

The UV behavior of this diagram is finite, as can be extracted from the list of propagators (185). Indeed, 
for large momenta, the corresponding integral of the diagram (1) behaves like ~ / A^pj^-^, which is 

perfectly finite in the UV. Therefore, limt^y (^^(x)w(>')) is not divergent at one loop. In the next sec- 
tion, we shall explicitly prove this. 

At two loops, it is not possible to present the same argument as there exists a diagram which can be 
logarithmically divergent: 



A 


95 











A ^ 



Figure 2: A possible divergent 2-loop diagram for (A^(x)(p9(y)). 
as can be checked from the list of propagators (185). 

-There is no quadratic coupling of CO and (5 to the gluon .sector, thus such a condensate would not directly influence the gluon 
propagator. 
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Secondly, we can also have a look at the mixing of the operators and (p(p. In the algebraic analysis, 
see appendix C, we have found that a mixing is possible between the different operators, see equation 
(173). This means that algebraically, a counterterm of the type QA^A^ is allowed. This counterterm is 
needed to cancel the infinities of the following type of diagrams: 




However, we can prove that there are no infinities at one loop, as the only possible diagram is given by, 




which is similar to the diagram in Figure 1. We can thus conclude that the mixing can only start at two 
loops. Again, we cannot exclude divergences at two loops, due to a similar diagram as in Figure 2. 



4 The effective action 

In this section, we shall try to calculate the effective action. The calculation is quite technical and shall 
therefore be split in different steps, although the result is reasonably compact and can be immediately 
found in expression (96). 

The energy functional can be written as 

g-w(G,x,G,G) ^ y[d4^][dc][dc][dZ;][d(p][d(p][dco][dco]e-^'^°^ (19) 

with EcGZ given by equation (18). We recall that in rf = 4 — e dimensions, we have the following 
dimensionalities, 

r n r , d—2 e 

_ ^^'^ _ ^ 
[g] - ~ 2 ' 

W = [el = [Gl=[G]=2, 

IQ = [a] = [x] = [p]=d-4 = -e. (20) 

4.1 The LCO formalism 

In order to calculate the effective action, we shall follow the local composite operator (LCO) formalism 
developed in [24, 25]. Let us outline the main idea. We start from a LCO O, in our case a local dimension 
two operator within a dimension four theory. As done several times, we couple the operator(s) of interest 
to an appropriate source(s) /, and add the term JO to the Lagrangian. This gives rise to a functional 
W{J) which we need to Legendre transform to find the effective potential. However, as already observed, 
novel infinities shall arise, which are proportional to J^. These infinities are due to the divergences in the 
correlator livcix^y {0{x)0{y)}, as explained in section 3.3. Therefore, in general, a term proportional to 
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is always needed in the counterterm, and the starting action needs to display a term-' ^i^. The novel 
parameter ^, called the LCO parameter, is needed to absorb the divergences in fi, i.e. 

5^y2 With the 

inclusion of the term the functional W{J) obeys the following homogeneous RGE 



(^^±+^{g^)^-^j{g^) I d*xJ^+n{g\Q^jw{J) = 0, (21) 
with T| , Q the running of ^, 

^.^^t; = ^{g\Q. (22) 

Notice that it is necessary to include the running of C, at this point. 

Now the question is, how can we determine this seemingly arbitrary parameter ^? This is possible 
by employing the renormalization group equations. We can write 

l^o4=^,-'{t^J^ + 8t^J^), (23) 

whereby the second term of the r.h.s. represents the counterterm. As the l.h.s. is independent from yu, we 
can derive both sides w.r.t. /u to find: 

-£K + 5Q+(^^,|-C+A-|;(5Q)-2Yy(g2)(C + 6Q=0, (24) 

whereby yj{g^) is the anomalous dimension of J. As we can consider ^ to be a function of g^, and by 
evoking the P function, 

^ig')=f^/ (25) 

the equation (24) becomes, 

p(g2)i^CU') = 2Y/(g')C+/(/) ■ (26) 
with /(g^) = e5^ — P(g'^) ^ (5Q + 2jG{g^)dt,. The general solution of this differential equation reads 

^(/) = + aexp 1^2^* |Mdzj , (27) 

with C,p{g^) a particular solution of (26). A possible particular solution is given by 

t^p{8^)^'^+cih + C2g^hU... . (28) 

whereby we have temporarily introduced the dependence on h. Notice therefore that the n-loop result 
for ^(p^) will require the (n+ 1) loop results of P(g'^), yj{g^) and f{g^). As we would like t, to be 
multiplicatively renormaUzable, we set a = 0. In this case we have that 



%')+5%') = Co=ZcC(/), (29) 

and we have removed the independent parameter a. Also, now that ^ is a function of g-^, the RGE (21) 
becomes 

(^M^^ + p(g2) ^ - Y,(g2) I d^x/^) W{J) = , (30) 



as deriving w.r.t. ^ is now incorporated in deriving w.r.t. g^. 

After determining the LCO parameter ^, the next step is to calculate the effective action by doing a 
Legendre transformation. However, it shall be easier to perform a Hubbard-Stratonovich transformation 
on W{J), whereby we introduce an auxiliary field o describing the composite operator O. In this way, 
we can get rid of the quadratic term in and a clear relation with the effective action emerges, as it will 
be shown later on in this section. We only need to mention that the case we are handling here is a bit 
more complicated due to the mixing of the operators Oi = (p,(p,- and O2 = AfjA^, and to the mixing of the 
vacuum divergences. However, the basic principles remain the same. 



'For an example, see the action (18), where the term — j t,x- — a.QQ — xQx is needed in the starting action. The sources Q and X 
are coupled to the LCO operators Oi = 9(9; and Ch = A^A^- Note that here, also a mixing term %Qi accounting for the divergences 
in Mm^^y {Oi{x)02,(y)) is present. 
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4.2 Differential equation for the LCO parameters ^, a, % and p 

We shall try to determine the four LCO parameters ^, a, % and p. We shall first derive a differential 
equation for these parameters, in an analogous way as in [24, 32]. As there can be mixing, we shall 
define 5^, 5(0 and 5% as follows 

1 , 



-Xa^-a^^Qi)-^^^Qa■^a■ 



while 6p can be defined independently: 

PoGqGo = it^ZgZ^pgg = ir^[\ 

We further define the anomalous dimension of G, 



pGG. 



^dn ^"^^ " ^^^^^^ 



i"g^G = -yG(/)G, 



(31) 



(32) 



(33) 



which is exactly the same as the anomalous dimension of G as Zg = Zg. To define the anomalous 
dimensions of Q and X, we start from equation (173): 

n 1 

(34) 







Zqq 




e 






Zt;Q Zx^ 




X 



a relation stemming from the algebraic renormalization. To the matrix Z, we can associate the anomalous 
dimension matrix T: 

fj^Z = Zr, (35) 

and thus 










0" 






Til 


Ytx. 



This matrix is then related to the anomalous dimension of the operators: 



Xo=ZX 



9Z„ ^ dX 
■0 = ti^X+Zp— 



dX 
% 



-rx. 



(36) 



(37) 



(38) 



so the anomalous dimensions of the sources Q and x is given by 

% 

With these definitions in mind, we can derive a differential equation for 5^, 5(n, 5% and 5p. We start with 
that of 5p. Starting from expression (32) and deriving w.r.t. /j, we find 



Q 




-IQQ ■ 




Q 


X 




.-r21 -Yxx 




X 



-e(p + 5p)+ ( n^p+M^{8p) ) -2yG(g')(p + 5p) = . 



(39) 



As we can consider p to be a function of g^, according to the standard LCO formalism, we can rewrite 
this equation as 



P(g2) |jp(/) = e(p + 5p) - p(/) A (5p) + 2yG(g2) (p + 5p) . 
As p is finite, we can even further simplify this into 

PU2)^p(g2) =2yG(g')p + e5p-p(g^)£.(5p) + 2YGU^)5p . 



(40) 



(41) 



In an analogous fashion, we can find the differential equations for 5^, 5co and 8%. If we derive (31) 
w.r.t. fj, we find the following set of coupled differential equations 



'dg^ 2 



|5^-^p(g2)|j(5Q+Yxx(g2)(^ + 6g, 



p(g2) ^a(/) = e6a - p(/) ^ (5a) + 2YeG (/) (« + 5a) + Fzi (/) (x + 8%) 



PU')9^5C(«') =£5x-P(g')3^(5x)+YeeU')(X + 5x)+YxxU')(X + 5x) + r2i(/)(C + 5Q • (42) 
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4.3 Determination of the LCO parameters 5^, 5a, 5% and 5p 

In order to determine the countertorm parameters 6^, 5a, 5x and Sp at one loop, we need to calculate 
the one loop divergence of the energy functional lV(2,x,G, G). The details of these calculations can be 
found in appendix E. From section 3.3, we know that at one loop, 5% should be zero. This observation 
shall serve as a check of our computations. 

In the appendix E, equation (202), we have found 

1 1 T 1 

Sex - -,4:?(A^^-l)^ 



5x = 

1 1 
£ Alt- 



5P = -l^iN^-if- (43) 



The value of 5^ provides already a first check of our results. In fact, this quantity has been calculated up 
to three loops, see [24, 33]. Our one loop value for 5^ coincides with that reported in [24, 33]. Secondly, 
we also see that indeed 5% = at one loop, which nicely confirms our diagrammatical power counting 
argimient. 

4.4 Solving the differential equations for ^, a, % and p 

In this section, we shall try to solve the differential equations (41) and (42), when possible. For these 
calculations, it is useful to keep in mind the (5 function, here given up to two loops 

P(g') = -e/-2(po/ + pig6 + o(/)) , (44) 



with 



in order to keep track of the orders. 
We start with (41), 

P(g2)^p(g2) = 2Ycte')p + e5p-p(/)g^(5p)+2Yc(g2)5p . (46) 

In order to solve this differential equation, we need to parameterize p as follows: 

p=^+pi+P2/ + 0(/). (47) 
g 

We also need the expUcit value of the anomalous dimension Jq. We have from the definition (33) that 



a 



yc{g')=p^\nZG, (48) 

and thus we need the value of Zq. From the renormaUzation factors (167) and (137), we find that 

YG(g') = -M^ lnZ<p = ln(Z-iz;'/2) . (49) 

In [34], the factors Zg and Z^i have been calculated up to three loops, 

7 - , I 13 1 iVg' , /-13 1 59 n f Ng^\^ 

^ ~ ^6el67t2^l 8 e2 + 16e; I 16;r2 j ' 



11 1 Ng'^ nil 1 _ 17 1\ / Ng 
~ 6" e 167? "^1,24^ £2 6"eyVl67t2 



2 \2 

(50) 
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So one can calculate yoig^) up to three loops if necessary. Here only the first loop shall be useful for our 
calculations, i.e. 



Yc(g') = 



3 Ng' 



+ . 



(51) 



4 16jt2 

as 5p, see equation (202), is only known up to lowest order. With this information, we can solve the 
differential equation (46) up to lowest order, by matching the corresponding orders in 



24 (N^-l) 
53 iVg2 



■+P1 + - 



(52) 



Unfortunately, we cannot solve the differential equation for pi as we would require the two loop value 
of 5p, which is however not easily computed. Therefore, in the current work, we leave this value as a 
parameter to be determined. 

Let us now turn to the set of differential equations (42). We can do a similar analysis as above for 
the first differential equation, namely 



2^ 9 



'3g2 2 2 
We shall again parameterize C, as follows: 



5C-^P(/)3^(5Q+Yxx(/)(; + 6Q 



;=^ + Ci+^2/ + 0(/). 
g 

In fact, we can even solve this differential equation to two loops. From [24, 33, 32], we know that 



(53) 



(54) 



5^ = 



163l2 



35 1 139 



1\ /g^N\ 



2 e2 6 e / V 16^2 / 



+ 



665 1 6629 1 

■ + - 



6 36 £2 



71551 231 
W + -16^(^) 



eJ\ 16ji2 J 



(55) 



and 



Zvz = l- 



35 1 



6 £ V 167l2 J 



2765 1 449 1 
^72"£2 ~ 48"£ 



+ 



V 16ji2 J 

113365 1 41579 1 



432 £3 ^ 576 £2 



+ 



75607 _^ 
2592 16^^ ' ' " 



\ 167l2 J 



SO that from (36) 



, 2^ 35 / g^N \ 449 / g^N \ ^ ( 



94363 9 

^^ + — a3) 

864 16^^ ' 



16jr2 



(56) 



(57) 



By solving the differential equation for ^, we can determine ^ to one loop order. In principle, we can even 
go one loop further with the known results. However, as we shall only determine the effective potential 
to one loop order, we do not need this next loop result. We find, 



A'2- 1 
167l2 



9 167t2 161 

+ 



13 g^N 52 



(58) 



see also [32]. 

The second and third differential equation of (42) are coupled. However, they can be simplified and 
decoupled as 8% = 0: 



p(/)^a(g^) = 2yee(/)a + £5a-p(g2)_(6a)+2yee(/)8(x + r2i(g^)X, 



P(g')3^5C(g') 



Yfie(/)x+YxT(g')3c+r2i(g')(C+5C) . 



(59) 
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Fortunately, we know that = at lowest order, from the diagrammatical argument in section 3.3. 
Therefore, we can set = + 0(g'*). When parameterizing as usual 

a=^+ai+a2/ + 0(/), X=^+li+lig' + 0{g^) , (60) 

we find for the solution of the differential equations 

2A{N^-lf 

= 35A^' 

Xo = 0. (61) 

4.5 Hubbard-Stratonovich transformations 

In this section, we shall get rid of the unwanted quadratic source dependence by the introduction of 
multiple Hubbard-Stratonovich (HS) fields. We can then rewrite the relevant part of the action in terms 
of finite fields and sources: 

t 

1^_„_. ^ „ 1 



d^x [zcZop - G(pf (p? + ZcZ^ - G(p>° + ZpZ2 pGG 



We shall now perform the following Hubbard-Stratonovich (HS) transformations by multiplying expres- 
sion (19) with the following unities^. 



1 



1 



y [doi]e" ^ i^^''{^-J+b^^"■^^-lXi^^^lH-t^r^l^d) 

J [d(J2]e" 4C'[4a'C'-x'^] ''"K ? +(6x'-2<.C')f ^/'A^-2cCV''^W+(4a'C'-z'^)^'-^/^e) ' 
/"r , , - 47k; /d''^(^ + j/^'"ZGZ^W+i*^/'ZGZ<pq)<p+Z|ZpP;U-^/^G+4ZpP,,-^AG)' 

/ [d03]e P gp ^ " ^ , 

/ [d04]e ^^f^o" ^ " ^ , (62) 

whereby we have introduced four new fields, Oi ,02. ^^3 and 04. By doing these HS transformations, we 
can remove the quadratic sources and rewrite the functional energy as 

g- w(ax,G,G) ^ I [j^^j [jp] [j^] [j^j ] [^i^j^] [daa] [do4] [dcp] [dip] [dco] [dm] 
with (|) = (A^,c,c,Z7,(p,(p,(B,ro) and 

+ ^ ^ + feX'-2aC e/2 ^^2 c e/2 

4^'[4a'^' - x'^] [4a'C' - X'^] 4a'C' - X'^ "^-^f 'P' 



''We dropped irrelevant normalization factors. 
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As these HS transformations do not put everything in the right form yet, we propose the following extra 
transformation 



2;' 2^' 



■>2 ■ 



(65) 



So (63) becomes 



X eL 



, (66) 



whereby 



?2 4a';'-x'2 g2 



2W-ax' £l„e/2^2 _ bx' -laC,' ^^2 _ ^X^„e/2 



+ 



4C'[4a'C'-x'2]' 



[4a'C'-x'2]^ W^') 4a'C-X'2 ^ ^' + 4ZpZ2 p + 



4ZpZGP g " 4ZpZcP ^ 

8 



Z2 
4Zpp' 



(67) 



Now acting with ^ | ^ and | ^ on the energy functional, before and after the HS transformation, 
gives us the following two relations, 

ZeeZ<p((pf(pf) + iz^z,;^(AX) = -n-'^'^^ , 



-ZaZxx(AX) 



e/2(Pl) 



(68) 



while acting with ^ _ and i _ 

°^ I G,G=0 I G,G=0 

ZcZp ((pep) 



or equivalently 



ZcZp ((pep) 



-£/2 (03+104) 
-E/2 - 104) 



(69) 



~e/2 (03) 
-E/2 (04) 



(70) 



4.6 The effective action 

If we introduce the parameters 

2 

m 

T 

P = 
nt - 



1 



4ocoCo-2xg 
1 

200^0 -Xo 

53iV 
48(A'2-1)2 

53Af 
48(A?2-1)2 



2 (2ao^(Ji-Xog02) 



2 (Xogoi - ^0^02) , 
(03 - i(J4)g , 



(71) 
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with (Xo, CoiXo given in equations (58)-(61), then the quadratical part of the Lagrangian (67) is given by 

- 4^ (? + ?) + :J,./>.2_MV'^W+.-iw.,-4«) . (72, 

We have left out the higher order terms as we shall only calculate the one loop effective potential . 

All details of the calculations of the effective potential have been collected in the appendix E. The 
final result for the effective potential is given by 

-3^.(^-.i.(^]-2(.^-l)|,.^3^(.^-l) 

1 48(iV^-l)V ,53 p, W 

2 53N \ ^ 24(]V2-l)2y ^2 



+ T3T2?-35Sr^^-^-l6^T-^«'+^""'- ^''^ 

whereby y\ , y2 andy3 are the solutions of the equation j"* + (m2 + 2M2)y2 _(_ _|_yj^4 _ ppt _^_2M^irP')y + 
m2A,4 + 1/2(p + ^ ^4^2 _ ^2ppt = and y4 and ys of the equation y'^ + 2M^y+M'^ - ppt = 0. 
We employed the MS scheme. 

4.7 Minimizing the effective potential to prove that the condensates are 
non-vanishing 

To simplify the calculations, let us set p = p^ = 0, which corresponds to the case of not considering 
the condensates (99) and ((p<p). For the moment, we are only considering (99), which already has the 
desired influence on the propagators, see the next section. In this case, the effective action simplifies, and 
becomes: 

gN^-ltn^ 24(N^-lfM'^ mN^-lm^ . ,,2 2 
+ T3T2?-35Sr^^-^l6^T-^«'+^""'- ^''^ 

whereby j2 and yj are are given by 4 ( —m^ — ± \/m'^ — 2M^irP- + M'* — 



In order to find the minimum, we should derive this action w.r.t. np- and and set the equations equal 
to zero. In addition, we should also impose the horizon condition (9). Therefore, we have the following 
three conditions, 

3r 3r 3r 

which have to be solved for M^, and . Unfortunately, it is impossible to solve these equations ex- 
actly due to the two unknown parameters ai and %\ . However, we would bice to know if the condensate 
(99) is present or not. For this, we need to uncover ifM2 = can be a solution of the above expression. 
We can strongly argue that this is not the case, and thus that / 0. 
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We shall start from expression (74) and derive w.r.t. M^, and X'^. As we would like to know if = 
can be a minimum of the potential, we further set = 0. We then obtain the following equations 



3 (hi (nP- - Vm4-4X4) - In (m^ + Vm^^^)) X'* 



43i2Vm4-4X4 
ll\/m''-4A,4(241n2- 17)m^ 



+ « Xi--T=0, 



Vm4-4^4 

+ 39 (-m'* + V'W*-4A,4m2 + 2X'^'J In Q (^m^ - v'm^^^) 
+ 39 (m'' + \/m'^-4XW - 2X'^^ In f ^ (^m^ + ^m'^-4X'^^ 



0, 



\/m4-4X4 
1 



9(VOT'*_4^4_^2jjjj/ 1 (^m^ - - 4^4 j j + 9 (^m^ + Vm* - 4:^,4^ 



ln(-(n? + \/m'^-4X'^U+Vm'^-4X'^{-l5+n6ln2) 



(76) 



whereby we have chosen to set^ j" = 2 and A' = 3. Now looking at the equation, we see that the second 
and third equation can be solved exactly for and X. There are even multiple solutions possible. We 
take the solution which has the lowest value for the effective action with M-^ = 0. However, for this so- 
lution to be also a solution of the first equation, these values should be very specific and the chance that 
they will also satisfy the first equation is practically non-existent, with a certain value of %i . Moreover, 
at a different scale p, the three equations will look slightly different. However, Xi is a number and stays 
the same. Therefore, it would be necessary that at all different scales these three equations can be solved 
exactly for only two parameters. This is practically impossible, leading to the conclusion that ^ 0. 
A similar reasoning can be worked out if p and/or p ' would be allowed. The main result is that it is 
impossible for all these condensates to be zero, making the associated refinement inevitable. 

hi conclusion, we have a firm indication that the condensate ((pep) is indeed present, thereby suggest- 
ing the dynamical transformation of the GZ framework into a refined GZ framework, with associated 
propagators that are in agreement with the most recent lattice data of [1, 2, 3, 4]. 

5 The gluon and the ghost propagator 
5.1 The gluon propagator 

The gluon propagator shall still be infrared suppressed and non-zero at zero momentum. Indeed, starting 
from the further refined action (14), the quadratic action is given by 



^quad 



2 f 



(77) 



whereby we have replaced the source x with mP, Q with -M^, g'^ with -6'jp and G'j with -S'^p^ and 
set all other sources equal to zero. From this, we can easily deduce the gluon propagator 



a;{p)aU-p) 



=. PuPv 

P 



2(M2+/)2_2ppt 



2M4/P-2//^-2M- (Ip"^ + -X-^ip + p') + 2in- ( (M- - p-) - - pp-) + 2p-(X^ - pp^') 



(78) 



' We work in units Ar^ = 1 . 
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with = Ig^N-f. If we assume that p = p^, we then find the following gluon propagator: 

M2 + p2 + p 



p'^+M^p^ + p^{p + m^)+m^ (m2 + p)+>.4 ' 



(79) 



which has exactly the same form as the refined gluon propagator (12). However, for the moment we 
cannot say whether p = p^ is the case or not. This shall be further investigated in [29] upon using lattice 
input. Notice that p, p^ as well as provide in an independent way that D{0) ^ 0. In principle, it could 
occur that M"* = ppt, giving D(0) = 0, but there is no obvious reason why this relation should have to 
hold. 



5.2 The ghost propagator 

The one loop ghost propagator is given by 



1 1 



with 



(27t)4 (k-q)^ 



(80) 



Ik^J (27t)4 



N^-lk^J (27t)4 {k-q)2 

2kfik\ f d'^q 1 



I 



k^ J {2%Y{k-qY 



r 



2(M2 + g2)"_2pp"'" 



2M4^2 + 2q(' + 2M2 (2^4 + ^4) - A,4(p + p') + Im'^ {{M^ + q^f - ppt) + 2^2(X' 



4 -ppt) 



As we are interested in the infrared behavior of this propagator, we expand the previous expression for 
small k^ 

d-1 f d'^q 1 



(231)'' 



2(M2 + ^2)2_2ppt 



2M4^2 + 2q6 + 2M2 (2^ + ^4) - A,4(p + pt) + 2m2 ((m2 + ^2)2 _ ppt) + 2^2(^4 _ ppt) 

+ 0(k^). (81) 

Let us now have a look at the gap equation. For this we can start from the (one-loop) effective action 
which can be written as (see the appendix E) 



(^-1)/ 



d^'q 
(2^ 



lnA + . 



with 



2M4^2_^296 + 2M2 (2?4_^^4) _:^4(p_^p1-)_^2m2 (^(m2 + ^2)2 _ ppt) +2^2(^4_ppt) 



2(M2 + ty2)2_2ppt 

and the . . . indicating parts independent from X. Setting A,4 = Ig^N^ , we rewrite the previous expression. 



f(l) 



Af2-1 d 



The gap equation is given by ^S- = 0, 



3^2 



d-l f d'^q 



2M2 + 2^2_p_pt 



2M4^2 + 2^6 + 2M2 (2^ + X4) - X4 (p + pt ) + 2m2 ( (M2 + ^2) 2 _ ppt^ _^ 2^2 (X4 - ppt) 



, (82) 
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where we have excluded the solution A, = 0. With the help of this gap equation, we can rewrite equation 
(81), 

d~l f d''q 



2M4^2 + + 2M2 (2^4 + X4) - X4(p + pt ) + 2m2 (m2 + ^2) 2 _ ppt^ ^ 2^2 (^4 _ ppt ) 

■^0{)?). (83) 

The integral in the above expression is finite. We can rewrite the integral as {d = 4) 



M^q^ +q<>+M^ {iq^ + ^4) - rX^ + nP-(^ (M2 + ^2 j 2 _ (^2 + ^2) j _,_ ^2 (^4 _ (^2 _,_ ^2)) 
with / = o(fe2 0) — 1, whereby we have parameterized 



p = r+is , 



p ' =r — is. 



(84) 



We further write 

/ = ^ ix (m^ - + s^) +x{M^ + r)) / (x^ +x^{2M^ + n?) 
Solution of cubic equation 

The next step would be to solve the cubic equation in the denominator of the equation above, 

+x^ (2m2 + m2)+x(M4 +2m2M2 + ^4 _ (^2 _^_j2)^_^^4(^2 _ +^2(^4 _ (^2 + ,2^) = . (86) 



In general, the roots are given by 



XI 



with 



X2 



X3 









a + 


T 








a + 


T 










a + 


T 



-1 +i\/3 ijm + y^ ^ —l — i^/3 ijm — 



-1 — i\/3 _3/w + -v/n — 1 +i\/3 a/m— -v/n 
^ 2^ 



2 - {n? -Nf) - 9 (,-2 + .v') j - 9 (,f? -M^ + 3r) , 

2(^m2_jy^2j ^^^2_^2j2_9p_^^2j^ _9^OT2_^2_^3^j;^4 

(;„2_^2^%3(^^^^2_^4y 



n2 



(87) 



(88) 



Of course, it is possible that two (or three) solutions coincide. This can be checked by calculating the 
discriminant 

A = ~4a^c + a^b'^ - 4b^ + ISabc - 27c^ . (89) 
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If A = 0, then the equation has three real roots and at least two are equal. 
Case 1: x\ ^X2^ x$ 

If 7^ X2 7^ X3, we can rewrite the integral / as. 



I = Ng' 



64jt2 



M'* - {r^ + s^) + {M^ + r)xi 1 



(jTl -X2)(xi -X3) x-xi 



i 



+ / dx 



M'^ -{r-^ + s^ ) + {M'^ + r)x2 



i 



+ / ck 



{X2-X3){x2-Xi){x-X2) Jo {X3-Xi){x3-Xi){x-X3) 
These integrals are now easy to solve, they all are of the type j dx^ = Inx. 



I = Ng' 



64jt2 



M'^-{r^ + s^) + iM^ + r)xi 
{xi-X2)(xi -X3) 



M'^ - [r- + S-) + (M- + r)x2 
(X2-X3)(x2-Xi) 



hl(x-Xi)|o 



(90) 



ln(x-X2)|o + 



{X3 -xi){x3 -xi) 



ln(x-X3)|j 



One could expect there is a problem at infinity, in contrast with what we have concluded before. However, 
as Ml + vi + wi = 0, the infinities cancel. We obtain, 



I = Ng' 



647t2 



M'^-{r^ + s^) + {M^ + r)xi 
(xi -X2)(xi -X3) 

M'^-{r^ + s'^) + {M^ + r)x2 



hi(-xi) 



, , , M4-(r2 + i2) + (M2 + r)x3, , , 
ix2-X3)ix2-X.) + (x3-.l)(x3-^t) ^"^-^^^ 



Case 2: xi = X2 7^ X3 

In this case, we can rewrite the integral / as 



-r 

Jo 
Jo 



647l2 



i: 



dx 



M^- (r2 + .v2) + (M2 + r)xi 1 



(xi -X3)2 



X — xl 



dx 



M'^-{r^ + s^) + {M^ + r)xi 1 



X-X3 



^ . ^ M'*-(r2+.s-2) + (M2 + r)x3 1. 

Jo XI -Xi {x-Xs)^ 



W2 



One can check that M2 + V2 = 0, so we can perform the integrations. 



I = Ng' 



647t2 



M"^ - (r- + s- ) + (M- + r)xi 



(xi -X3)2 

M4-(r2 + i2) + (M2 + r)xi 
(xi -X3)2 



In(-xl) 



ln(-X3) - 



M'^-{r^ + s^) + {M^ + r)x3 1 



X1-X3 



Xj 



(91) 



(92) 



(93) 
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Case 3: xi =X2= x^ 



Finally, in this case we can write 
so after integration 



(M^ + r) rdx-. — ^--^ + {M'^-{r^+s^) + {M^ + r)xi) H dx-. — ^- 
Jo (x-xiy Jo (x-: 



I = Ng' 



64ji2 



m2 + r _^ Af4 - (r2 +s2) + (M^ + r)xi 



XX 



lx\ 



(94) 



(95) 



Now we can make some conclusions. Looking at the different cases, it looks almost certain that / 7^ 0, as 
very specific values of the condensates would be needed to take care of this. Therefore, we have strong 
indications that the ghost propagator is not enhanced, in addition to the nonvanishing gluon propagator 
at zero momentum. 



6 Conclusion 

Although this paper is quite technical, the conclusions are quite simple. Firstly, we have shown that using 
the GZ action, more condensates can influence the dynamics. We have investigated in detail the following 
condensates: (A]^°) , (cpf 9° ) , (cpftpf) and ((pf (pf ) whereby the latter two were never investigated before. 
We have proven that we can renormaUze the GZ action in the presence of these condensates. In particu- 
lar, a renormaUzable effective potential, compatible with the renormaUzation group, can be constructed 
for the associated local composite operators. 

Secondly, for the first time, we were able to calculate the one loop effective potential in the LCO formal- 
ism: 

-2(M2 + ppt)]+^[-^(.^-2.V.V?ln(f)+.lln(^+.iln(^ 

1 48(A-^-l)2 / 53 pi W 

^ 2 53iV ^ 24(^2-1)2; g2 

9A^2_i 4 24 (7V2_i)2^4 161Af2_lw4 

whereby y 1 , yi and y-x, are the solutions of the equation >-3 + (m2 + 2M2)>.2 + (X" - ppt +2M2m2)y + 
hf-X^ + l/2(p + p"'' ) A,"* + M^n? - m2ppt ^ and yi, and y^ of the equation -|- 2M2y -i-M'' - pp''' = 0. 
Unfortunately, due to the existence of yet unknown higher loop parameters, i.e. tti, pi and Xi, in the one 
loop effective action, we are yet unable to give an estimate for the different condensates. Nevertheless, 
we have been able to already provide strong indications that some condensates are in fact non-zero and 
shall lower the effective action. We hope to come back to the explicit computation of the parameters ttj, 
Pl and Xi in the future. In particular, one should compute the divergences of the vacuum diagram in Fig- 
ure 2, the similar one for the mixing, and other divergent 2 loop diagrams steimning from the operators 
(pep and 99. Once this task will be executed, all information is available to actually work out the one loop 
effective potential and to investigate its structure and the associated formation of the RGZ condensates. 

Thirdly, we have also shown that in this further refined framework, the gluon propagator is non zero 
at zero momentum, and the ghost propagator will be non-enhanced. 

A complementary approach to the current one, is to find out to what extent a gluon propagator of the 
type (78) or ghost propagator of the type (80) could describe the lattice data, not only qualitatively, but 
also quantitatively. This is current under investigation in [29] for different space time dimensions. In [4] 
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it was already shown that a RGZ propagator (78) reproduces the SU(3) data very well. 

Another question which was not answered here, is whether a(fc^), see equation (80), is in fact smaller 
than one. This is necessary in order to be assured to stay within the Gribov horizon. However, this 
question shall also be addressed in [29], and we refer to this paper for further details on this matter 
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Sy = -y'gf d'x (f''A"y;+r'"A"^t; + ^{n^~ i) y") ■ (98) 



A Recapitulation of the Gribov-Zwanziger action and of its 
renormalizability 

In this appendix, we shall repeat the complete proof of the renormaUzation of the Gribov-Zwanziger 
action [28]. 

A.l The Gribov-Zwanziger action and the BRST symmetry 

We start with the Gribov-Zwanziger action, 

Sgz = So + Sy, (97) 

with 

So = SYM + Sgi + I d''x (yfd^ (of <pf ) - mfd^ (Df4) - gf'd^mfD^/Ci^) , 

d 

8 

We recall that we have simplified the notation of the additional fields ^cpj]'^, cp^^'.mJJ'^, co^'^j in So as So 
displays a symmetry with respect to the composite index ; = Therefore, we have set 

((p,r,(p-.05-,a)-) = ((pf,(pf,5)f,cof) . (99) 

The BRST variations of all the fields are given by 

■s*" = 0, 
jcof = 0, 

s^f = 0. (100) 

However, due to the y dependent term, Sy, the Gribov-Zwanziger action breaks the BRST symmetry 
softly [20, 22], see eq. (1 1). In order to discuss the renormaUzabiUty of Sqz, we should treat the breaking 
as a composite operator to be introduced into the action by means of a suitable set of external sources. 
This procedure can be done in a BRST invariant way, by embedding Sqz into a larger action, namely 

I<3Z = Sym + % + 5o+Ss, (101) 

whereby 

(102) 

We have introduced 3 new doublets ([/«', M°'), (V/, Np and (T^', R^') with the following BRST 
transformations, and 

sT^'=R"J, sl^' = 0. (103) 

We have therefore restored the broken BRST at the expense of introducing new sources. However, we 
do not want to alter our original theory (97). Therefore, at the end, we have to set the sources equal to 
the following values: 

Trail ^ ^ jail ^ q 

Iphys Iphys Iphys 

<L =-<L ='/'5"'V- (104) 

^ Iphys Iphys ^ Iphys 















5®f 


= 9?, 
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Table 1: Quantum numbers of the fields. 







c" 


c" 


b" 






(Of 


ml 


dimension 


1 





2 


2 


1 


1 


1 


1 


ghost number 





1 


-1 











1 


-1 


2/-charge 














1 


-1 


1 


-1 



Table 2: Quantum numbers of the sources. 





jural 


V 


yai 






K- 


L" 


2 


2 


2 


2 


2 


2 


3 


4 


-1 





1 








-1 


-1 


-2 


-1 


-1 


1 


1 


1 


1 









A.2 The Ward identities 

Following the procedure of the algebraic renormalization outlined in [35], we should try to find as many 
Ward identities as possible. Before doing this, in order to be able to write the Slavnov-Taylor identity, 
we first have to couple all nonhnear BRST transformations to a new source. Looking at (100), we see 
that only A", and c" transform nonlinearly under the BRST s. Therefore, we add the following term to 
the action Lqz, 

Sen = ld''x(^-K"^{D^cy+^-gL''f>''c''c'y (105) 

with and L" two new sources which shall be put to zero at the end, 

^;lphys = ^''lphys = 0- (106) 
These sources are invariant under the BRST transformation, 

sK^ = 0, iL" = . (107) 

The new action is therefore given by 

Sqz = ^Z + 5ext. (108) 

The next step is now to find the Ward identities obeyed by the action Lq^- We have enhsted all the 
identities below: 

1. The Slavnov-Taylor identity is given by 



5(E^}z)=0, (109) 



with 



2. The U (/) invariance is given by 

[/yE^z = 0, (110) 



By means of the diagonal operator Qf = Uu, the /-valued fields and sources can be assigned an 
additional charge. One can find all quantum numbers in Table 1 and Table 2. 
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3. The Landau gauge condition reads 



4. The antighost equation yields 



5c« " 8Ki 
5. The linearly broken local constraints yield 



^ = 3,A^. (Ill) 



^''o^+a,^ = 0- (112) 



jZ, I -\ oz I „^ Tui uz rt-^at'CyiOi/a 

GZ I ~\ GZ ^jrahcjr-h GZ _ jrahc Ah j jci 



6. The exact %,] symmetry reads 

%-2gz = 0, (114) 

with 

7. The integrated Ward identity is given by 

/"d'^xl c^^^+ro^^^+i/"'^^ 1 =0 (116) 
J y 5(0? + ' 5c« + " 5Z« j ^ ^ 

Here we should add that due to the presence of the sources 7^' and the powerful ghost Ward identity 
[35] is broken, and we are unable to restore this identity. For the standard Yang-Mills theory, this identity 
has the following form 

^"(5ym + %)=A^i, (117) 



with 



and 



A«i = ^1 d'xf' (kX - L'c') , (119) 

i.e. a linear breaking. However, it shall turn out that this is not a problem for the renormalization proce- 
dure being undertaken, see later. 



A.3 The counterterm 

The next step in the algebraic renormaUzation is to translate all these symmetries, which are not anoma- 
lous, into constraints on the counterterm Eq^' which is an integrated polynomial in the fields and sources 
of dimension four and with ghost number zero. The classical action Eq^ changes under quantum correc- 
tions according to 

Sgz->S^z + ^Sgz, (120) 

whereby h is the perturbation parameter. Demanding that the perturbed action (Lq^ + '^^gz) fulfills the 
same set of Ward identities obeyed by IIq^,, it follows that the counterterm is constrained by the 
following identities: 
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1. The linearized Slavnov-Taylor identity yields 

«£&z = 0, (121) 
with !S the nilpotent linearized Slavnov-Taylor operator. 



and 

2. The U (/) invariance gives 

3. The Landau gauge condition 

4. The antighost equation 



+ ^^5W^"^5^+<6Zf +<5^+^^'675?)' ^^^^^ 
= . (123) 

tZ/yEoz = 0. (124) 



8b" 



(125) 



^+3.5f = 0. (126) 



5. The linearly broken local constraints yield 



8 -,8 abc—h 8 



S.^-^r'sf^l^GZ = 0. (127) 



6. The exact symmetry imposes 



%2gz = , (128) 



with given in (1 15). 
7. Finally, the integrated Ward identity becomes 



The most general counterterm Zq^ of d = 4, which obeys the linearized Slavnov-Taylor identity, has 
ghost number zero, and vanishing Qf number, can be written as 



£gz = aoSYM + ^ J d^'x|fliA-;A^ + fl29;.c"A^ + fl3rc«+a4(/;'9^(pf + a5V;"aX 

+hgfabcRt^'lAl + hgfabcT^'tiAl + ^'5<9;,©f + ^e^^'a^^Cpf } , (130) 

with flo, • • ■ jfll 1 5 • • • arbitrary parameters. Now we can impose the constraints on the counterterm. 
Firstly, although the ghost Ward identity (117) is broken, we know that this is not so in the standard 
Yang-Mills case. Therefore, we can already set 03 = as this term is not allowed in the counterterm of 
the standard Yang-Mills action, which is a special case of the action we are studying^. Secondly, due to 
the Landau gauge condition (3.) and the antighost equation (4.) we find, 

ai = 02 ■ (131) 

particular, since we will always assume the use of a mass independent renormalization scheme, we may compute with all 
external mass scales (= sources) equal to zero. Said otherwise, 03 is completely determined by the dynamics of the original Yang-Mills 
action, in which case it is known to vanish to all orders. 
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Next, the linearly broken constraints (5.) give the following relations 

ai = -ag = -09 = «io = an = -bs = bi , 

04 = 05 = — flg = fl? ; bi = b2 = b^ = b^ = . (132) 

The Rij symmetry (6.) does not give any new information, while the integrated Ward identity (7.) relates 
the two previous strings of parameters: 

ai = —as = ^"9 = '3!io = an = — ^3 = = 03 = 04 = —05 = ag . (133) 
Taking all this information together, we obtain the following counterterm 



A.4 The renormalization factors 

As a final step, we have to show that the counterterm (134) can be reabsorbed by means of a multiplicative 
renormaUzation of the fields and sources. If we try to absorb the counterterm into the original action, we 
easily find. 



V 1 1,^0 

Zg = 1— n— . 



(135) 



and 



zk = zy^ 

Zl = (136) 

The results (135) are already known from the renormalization of the original Yang-Mills action in the 
Landau gauge [35]. Further, we also obtain 

-1/2 -1/4 _ 01 



-1/2-1/4 
^A ' 



z-1 

-1/2-1/4 



■ Zg ^/^Z^ ^^'^ . (137) 

This concludes the proof of the renormaUzabiUty of the action (97) which is the physical limit of Zq^- 
Notice that in the physical limit (104), we have that 

Zf = Z'^l'^Zl^l'^ . (138) 



71/2 


= 4''=z. 


71/2 


7-1/2 

— ^A ' 


7I/2 


- z-1 


Zm 


_ J "^1 _ 

2 


Zn 


7-1/2 

— ^A ' 


Zu 




Zv 


1 J, ''I 

= i-''T = 


Zt 




Zr 


= l-h"-^- 
2 
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B Inclusion of the operator in the Gribov-Zwanziger ac- 
tion 

For the benefit of the reader, let us also repeat the renomialization of the operator in the Gribov- 
Zwanziger action, which was first tackled in [36]. In this paper, it was shown that the presence of the 
condensate (A^) does not spoil the renormalizability of the GZ action. The GZ action with inclusion of 
the local composite operator AJ^A" is given by 

■Sagz = 5'gz + 542 , (139) 

whereby 

= ' (140) 

with 1 a new source invariant under the BRST transformation .v and ^ a new parameter. The renormaliza- 
tion can be done very easily with the help of the previous section. 



B.l The starting action and the BRST 

Again, we shall make 5agz BRST invariant. We define 

^AGZ = ^GZ + ^A2 

whereby TIqj. given in expression (108) and 



-xA^«+llA2a,c''--?T2 



(141) 



(142) 



with T] a new source and .sT) = x, so that (r| , t) forms a doublet. At the end, we replace all the sources 
with their physical values, see expression (104) and (106), and in addition 



■Plphys '■ 



= 0, 



(143) 



so one recovers 5agz again. 



B.2 The Ward identities 

It is now easily checked that the Ward identities 1-7 of section A.2 remain preserved. Obviously, the 
Slavnov-Taylor identity receives an extra term, 

5(I:agz)=0, (144) 

whereby 

jd I SSagz 5£agz . SEagz SEagz 



r^^v \ f Ad ( oSaGZ O^AGZ , oZaGZ oZaC 
^(SaGz) = Jdx y-J^^ + ^^^a 



^'^ 5c« mf 6(pf + " 8y«' + " 8r/ +^ 5ti 

B.3 The counterterm 

As all the Ward identities remain the same, it is easy to check that the counterterm is given by 

£agz = i:^,z + /d^^(f^A^2 + y^^' + («2-fli)TlA23,c'') . (145) 

whereby counterterm (134). This counterterm can be absorbed in the original action, Eagz 

leading to the same renormaUzation factors as in equations (135)-(137). 
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1/2 

In addition is related to Zg and Z^' [36] 



Zj — ZgZ 



-1/2 



and and Zq are given by 



1 + h{—aj — 2a2 + 4ai — 2ao) . 



(146) 



(147) 



C Renormalization of the further refined action 



C.l The starting action 

Let us repeat the starting action (14), 

^CGZ = ^GZ + + ^<l><t + ^ma + "^fipjSof + '5([Kp,0)(p + '^vac i 
whereby Zq^ is given by equation (108), Z^2 by (142) and 

J d!^xs{Vmf(of) = J dt^xlwrnfaf -V^fcof] , 
^ J d'^xs{G'%f^"j) = J d-^x 

J d'-x \k{G'^W - IH'W) + - IWV''^ 

-J d'^x[a{QQ + QW) + P{QW + WW)+xQt: + 8Wt:] . 



J(pq),oxp 



-G'J(?1(?''j-H'j&f((>''j 



(148) 



(149) 



C.2 The Ward identities 

With the help of appendix A, we can easily summarize all Ward identities obeyed by the action Scgz 
1. The Slavnov-Taylor identity reads 

5(IcGz)=0, (150) 

with 



^(EcGz) 



f 4 f 5ZcG 



SScGZ S^cGZ , SZcGZ SECGZ , ,a 5^GZ , -a S^CGZ 



8A« 



5L" 5c« 



smf 



co?^ +M--fe +iV--fe + Q^^°^ 



^ 57^' 



5y 5ti ^q'} dH'J J 



2. For the U{f) invariance we now have 



whereby 



0, 



(151) 



r/.._ f d4v((D?_^_(D«— +m?— -m"— -M"^'— ^ Ifj ^ I Af«' ^ 

U^j-JdX ^(p, (p, + CO, (0, g^^. (7^ g^^^. +7V^ . 



^^^^^^^ 



-2G 



, 6 



By means of the diagonal operator Qf = Uu, the single j- valued fields and sources still turn out to 
possess an additional quantum number. 
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3. The Landau gauge condition and the antighost equation are given by 



;gz 



Ml > (152) 



-^^^^^ = ^'^'^ 

4. The linearly broken local constraints yield 

5ScGZ 8ZcGZ , W/STcGZ _ „rabcj,b^,ci, 

^^^"^-^■^ '^'^b^ - (154) 

whereby the ... are extra linear breaking terms irrelevant for our purposes. 

5. The exact l^ij symmetry is broken beyond simple repair. 

6. The integrated Ward Identity is broken also beyond simple repair. 

7. There is however a new identity: 



§ScGZ _ 8ScGZ 
5P ~ SV 



(155) 





AO 


c" 


c" 


b" 




9^' 




col' 


U"' 


^^^^ 






dimension 


1 





2 


2 


1 


1 


1 


1 


2 


2 


2 


2 


ghost number 





1 


-1 











1 


-1 


-1 





1 





^/-charge 














1 


-1 


1 


-1 


-1 


-1 


1 


1 











U' 


G 


p 


w 




T 


11 


G'' 


G" 


//'■'■ 




dimension 


2 


2 


3 


4 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


ghost number 





-1 


-1 


-2 





-1 





-1 





-1 








-1 


1 


Qf -charge 


1 


1 


























-2 


2 


-2 


2 



Table 3: Quantum numbers of the fields and sources. 



C.3 The counterterm 

These identities (150)-(155) can be translated into constraints on the counterterm according to the quan- 
tum action principe (QAP), see [35]. Unfortunately, many identities are broken due to the introduction of 
these d = 2 operators. However, we are using mass independent renormalization schemes and therefore, 
the new massive sources {P, Q, V, W , G'^ , G'"', H'j, h'^) cannot influence the counterterm of the original 
GZ action (134) since they are coupled to <i = 2 operators. Said otherwise, there are no new vertices 
capable of destroying the UV-structure of the original GZ theory (134). We only need to check whether 
these operators themselves are renormalizable. Thus, the counterterm is given by 

^CGZ = ^GZ + + ^P-H 1 (156) 

with given by equation (134), and ZJ^ given by 

^1 = 1 d*x (yXA^2 + 1^x2 + (fl2 - aOiiA^a^c") , (157) 

as already determined in (145). IFp ^ is dependent of all the sources (P, Q, V, W, G'^, W, H'J, H'^), 
is of dimension 4, ghost number —1 and 2/ = and obeys the remaining Ward identities. Due to the 
Unearly broken constraints we find 

-^=^=0, -^^=0, -^^=0, -^^=0. (158) 

3(p 3(p 3co 9(0 
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Therefore, 



I.'p_n = 'Sz J d^x {biPA"/L"^ + b2VA"/L"^ + biQP + b4QV + bsWP + heWV + bjP'c + biiVT: 

+ bqQn + ftioWTi + ciH'jW + C2H''G") , (159) 
whereby bi, . . ., C2 are arbitrary constants. By invoking the new identity 

we can write 

Z^_^ = J d*xbi[{Q + W)A"^"^ + 2{P + V)df,c"A"^]+b3QQ + b4QW + b6WW + bjQx + bfiW': 

+ ci{G'jG'^ -2H'jH'^)+C2{G"G-'-' -2H"W^) . (161) 
Let us notice that due to the U (/) constraint, the term in C2 is only present when 

QUa^q ^ 2HPPW = GifG'' + IWV^ , (162) 
which is indeed the case due to hermiticity. 

C.4 The renormalization factors 

Let us now try to reabsorb this counterterm into the starting action (14). We shall split this analysis into 
three parts, according to 

+ = IL'j + ZJ, + , (163) 

whereby 

Zj = J d'^x ci {G'^G'^ - 2H'V') + C2(G"G'^ - 2H''H") , 
Zjj = |d4xi>i[(e + W)A^°+2(i>+V)9^c''A«] + |TA^« + («2-ai)TlA^a^c", 
IFjii = Jd^xb3QQ + b4QW + b6WW + bjQi + bsWx+"^C,x^ , (164) 
are the three parts which we shall try to absorb separately. 



Firstly, we start with the vacuum counterterm connected to the arbitrary parameters ci and C2. If we 
redefine ci and C2, we can write 

Zj = j d'^x ciK{G'jG'-' - 2H''W) + C2KG"G'^ - 2H"H") , (165) 

and if we define 

H'^=ZjiH'', H'Q^=ZHH'j, G'1=ZjP'' , G^=ZgG'^, Ko=ZkK, Xo=Zx^, (166) 
we find for the renormalization factors of the new sources and the LCO parameters k and X: 

7_ _ 7-1/27-1/2 

Zk = (l + ci)Z^lZc' = (l+ci)Z^lz-l , 

Zx = (l + C2)Z=iZei = (l+C2)Z^iz^i , (167) 
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and thus the part Hj can absorbed in the starting action. 
Secondly, let us focus on Ljj 

Z'u = I d-^x bi HQ + W)AlAl + 2{P + V)d^c''A''^ + y xA^^ + (az - ai ) ^A^^dj.c" . (168) 
We propose the following mixing matrix: 

/ 2qq Zqw Zqi 





ZwQ Zww Zw^ \ \ W \ . (169) 

\ z^e z^w z^x 

• From 

eo<,o9?.o = \2qqQ + Zg^iy +Ze,x]z^(pf (p? = q<s}1^i , (i70) 

we find that Zqq = Z^ \ while Zqw = Zq^ = 0. 

• From 

WoS^oKo = [ZwqQ+ZwwW +Zwxx\Z^^(sfi = lV(p?(pf , (171) 

we find that Zww —^^< while Zwq = Zwz = 0. 

• Finally, from 



l^oAl^fAlfi = \[ZxqQ + Z^wW +Z,x1]ZaAIAI 



2 

we obtain Zxx = Zx = (1 +a2)Z^\ andZ^g = Z^w = 2^1 
In summary, we find the following matrix 



^ {\+a2)xA''^l + biQA''^l + bxWAlAl , (172) 



Go \ 












Wo 




ft 







to / 




V Zt:w 


ZxW 


Z-tt 




(173) 



Now that we have the mixing matrix at our disposal, we can pass to the corresponding bare operators by 
taking the inverse of this matrix. 



(174) 




2(p 







( Qo 







1) 


Wq 








\ to 









Subsequently, we can derive the corresponding mixing matrix for the operators, since insertions of an 
operator correspond to derivatives w.r.t. to the appropriate source of the generating functional Z'^{Q,W,x). 
In particular. 



1 2 _ 8Z'{Q,W,x) 



8x0 



Xo=0 



86SZ^(e,W,x) ^ m 6Z'{Q,W,x) ^ 5x 5Z"(2,iy,x) 



5xo 5g 5xo SW 5xo 5x 
^Al = ^a\ (175) 

and similarly for (pfoVi'o ^"'^ '*^;'o'^?o- '^^^^ ^^^^ t° transpose of the previous matrix. 



W,o9-o 

'^0 





(176) 



We can make some observations from this matrix. Firstly, we find that Aq does not contain the operators 
(p"(p" and ro"co°. This is already a first check on our results as without these latter two operators the 
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GZ action including is renormalizable, as we have shown akeady in the appendix B. Secondly, we 
observe that 



(177) 



meaning that the mixing with A disappears again when recombining the two operators in a certain way. 
In fact, this is the operator (^f (pf — ro^co^') which we have investigated using the RGZ action [22] and no 
mixing with appears for this operator. 

We can do a completely analogous reasoning for the part in '^^l<fA^. We first set 7 + P = X. We propose 



From 



110 



-(^o)[<P^:o<.o] 



Zqx 



(178) 



-[ZxxX+Zx^Tl]zi/'z^/'(p?< 



we find that Zxx = Z^ ^'^^Zoj while Zxti = 0. 
Also, from 

-^^U^A = \Z^X + Zp^Ti]z]/^zy^A^a^c" = ( 1 + «2 - ai ) f\^l\c'' + 2fciXA;29^c° . 
we obtain Z,,ri = Zq = ( 1 + ai — oi ) Z^ '^'^Zc and ZpA- = 2Z)i . 



Therefore, we find that 



1/2 1/2 
Z/ Zc' -2/71 







Zn 



(179) 



Again, we find that q^/^cq does not contain (pJ'qCO^q, which is necessary as the GZ action with the 
inclusion of is renormalizable. We also see that, when setting V = —P, X = 0, the mixing with A^ 
disappears again. 



Thirdly, the vacuum term EJ^^ has the following form 



(180) 



we know that setting Q = —W has to return the vacuum term from the RGZ action ~ 0421+ ^t,x^. 
Therefore, we may set 



^3-^4 + ^6 = 0. 



(181) 



In this case, the vacuum term reduces to 



-cia(QQ^QW)-c2^{QW + WW) - c^xQx - c^mx + -j- ^x^ 



(182) 



where we have extracted a, P, % and 8 and some minus signs for convenience. If we allow mixing 
between the different parameters. 



( 

Po 
Xo 
80 
V Co / 



/ Z„|3 Z^i Z„g Z„^ \ 



pa ^pp ^Px ^PS 



^Z5 ^X^ 



Zga Zsj, Zgg Zg^ 



( a\ 
P 
X 
5 

V C / 



\ Z^p / 

when absorbing the counterterm, we find for the mixing matrix of the LCO parameters 





Zap 


Zax 


Za8 




Zpa 


Zpp 


Zpx 


Zp8 


zp^ 


^xa 


ZxP 


Zxx 


Zj,g 


ZxC 




Zap 


Zsx 


ZgS 


Zg^ 


V 


z^p 


h% 


Z^8 







l+ci 















1+C2 





ZggZtiv 






Zqq 










1+C3 
















1+C4 

ZqqZzx 


v 















2Z' 



00 



ZrwZ: 



^^QQ 
ZrH/7^; 

_ 7xH"^;; 

, ^ee 

1— fl3 



In summary, we have proven the action (148) to be renormalizable. 



(183) 



(184) 
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D List of propagators 



We give here the Ust of propagators which can be calculated from the GZ action (2): 



with 



Al{p)bHk) 



{b''{p)b\k) 



5"<^5*''5^^6(p + yk) (231)4. 
8«''^5(p + fc)(23i)4, 



r,4_ 



■ 1_ 

' ^2 



-i^?f''h{p + k){2nf , 



P 



5«*_5(p + *:)(23i)\ 



X. 

p' 

A",{p)%{k)) = r''^-^^P^{p)i2n)'8{p + k). 
b"{p)%{k) 



fbc;p^^(2Tlfb{p + k). 



1, 



8'"'8'"'8^]{2nyd{p + k)., 



% {P)% {k) 



fubr rcdr n 



p^ip'^ + lg^Nf) 



{2nyd{p + k) . (185) 



(186) 



E Details of tlie calculation of the effective action for the fur- 
ther refined GZ action 

E.l Determination of the LCO parameters 5^, 5a, 5% and 5p 

We shall start from expression (18), determine the quadratic part, and integrate out all the fields. The 
quadratic action is given by 



I d\ [ecpf 9f + I'^A^Al - i ^x^ - aQQ - xgx] + / A'x 



whereby we have immediately integrated out the ghost fields, c,c, (jO,c5, as they only appear trivially. We 
have also already integrated out the i)-field whereby a is formally equal to zero. 

As a first step, we integrate out the cp and ^ fields. For this, we shall spUt <p, ^, G and G into real 
and imaginary components: 



G = X + iY , G = X- \Y , 

so that the part depending on (p and (p in expression (187) becomes 



(187) 



-2YU^V^ + pX^ + 



2{d^ + Q + X) -2Y 

-2Y 2{d^ + Q-X) 



yah 



be • 
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Therefore, applying Gaussian integration, we find for the integration over <p and (p 



|[d(p][df]exp[-Z*g^]=exp 



2^ " 



d^ + Q-X 
84 + 2e92 + e2_x2-y2 



(detP;^'''^'')-V2 



whereby we recall that X is defined as A,'* = 2"fg^N. is given by 



'ftV>- 



2{d^ + Q + X) 
-2Y 



-27 
2(32 + 2-^) 



(188) 



(189) 



and the . . . stand for the other terms in S^GZ' (l^^)' i-®- terms purely in A and the vacuum terms. The 
second step is to integrate out the gluon field Aj^. Combining the expression (188) with the terms purely 
in A from the quadratic action, we obtain. 



det ^-6^3- + ^1 - a J - 94 + 2532 + ^2 _ ^,2 _ ^2 J + 



-t5. 



-1-1/2 



(190) 



Therefore, the total effective action at one loop is given by 

^-WiQ,.,Gfi) ^ (det<-'')-l/2 [det(-5,v92 + (1 - a) ^^^^ 



a4+2ea2+G2_x2-y2 



-1/2 



Si^x\- 1 ^T^-aee-xgrl +pGG]] 



(191) 



In order to find 5^, 5a, 5x and 5p at one loop, we need to find the first order infinities of the previous 
expression. These shall be present in the two determinants which we need to evaluate. 

Let us start with the first determinant of F^''^^. In general, we can write 



(det<'^'')-i/2=e-2 



^g_id(w2_l)2TrlnP (192) 



As we are taking the trace, we know that TrlnP = TrlnP' with P' the diagonalization of P. Therefore, 
after diagonalization, we find 



{detP^'"')-'/^ 

= exp -^d{N^ - l)^Tr (ln{-d^ -Q+Vx^ + Y^)+ln{-d^ -Q- Vx^ + Y^fj 
Employing the standard formula, [37] 

Trln(-a2+M2; 
we obtain the following infinity 

'1 (A'2- 1)2 



n-d/2) 1 

(47t)''/2 (M2)-''/2 



(deti')-V2 
whereby ci is a constant term. 



exp 



43l2 



[e2+x2 + y2j +ci 



(193) 



(194) 



(195) 



The second determinant requires a bit more effort to be evaluated. Let us call the corresponding ma- 
trix K. We thus calculate 



(det^g 



ah\-l/2 . 



e 2 



j(iV^-l)TrlnJr^ 



(196) 



Therefore, we need to determine 



TrlnA^uv = Trln | 8, 



+ x 



Trln 6, 



a4 + 2G92 + e2_x2_y2 
1 

( 92 x*( I A 



(197) 
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For the first term, we can easily take the trace over the Lorentz indices, while for the second term, we 

.2 

need to use ln{l + x) = x — 'j- + ■ ■ then take the trace of the diagonal elements of the second term, and 
again employ x— ^ + . . . = ln{l +x). After these operations, we obtain 



TrlnA"^ = rfTrln ^ ^- 



which can be written as 



d'^ + Q-X 



+ Trln f 1 + -. \ r ^ f 1-- 

V y O A, \^34+2|2a2+g2_x2_y2 



muK„ . w- i)Tri„ ((-3' - X' ( e.^2ey+g-r-y^ 

The first term of this expression can be written as^ 

(^^ - 1 ) [Trln(p6 + _ 22) + ^X'^ + g2 _ ^2 _ y2 _ ^q^^ p2 _ q^A ^ ^yA ^ q2^ 
- X\ - Y\) - Trln - 22^2 + g2 _ ^2 _ y2^ j 
=[d-\) fTrln(p2 +Trln(p2 - xt) +lv\n{p^ - x^) -IxXnip^ - x^) -TY\n(p^ - xs) ] , (198) 



whereby x\ , X2 and x-i are the solutions of the equation + (x — 22)^2 + {X^ + —X^ — Y^ — 2Qx)x — 
QX'^ + XX^ + Q\-X\- Y\ = and X4 and X5 of the equation ^2 - 22x + 2^ - A"^ - 7^ = 0. After 
determining xi , . . . ,X5, we can apply the standard formula (194) again, so we ultimately find for the first 
term 

3 1 



163i2 e 



(x2_2X4)+C2, (199) 



with C2 a constant, which is not of our current interest. For the second term of (197), we can perform an 

analogous analysis, whereby we find that this term is proportional to a and therefore does not contribute 
to the determinant as a — > 0. Therefore, the second determinant ultimately gives: 



(detZ^)-i/2 = exp 
We can now combine both results (195) and (200) to find 

(yv2-l) 1 /3 , , o ,,,^2,^7^x 3 



(200) 



W{Q,x,G,G) = -^-^^-{^-i^ + {N^-l){Q' + GG)-'-X'\+c, (201) 



with c a constant term. Therefore, at one loop we obtain 

1 1 o 

5« - -e4^^(^^-l)^ 



5x = 0, 

5P = e4^(^'-l)'- (202) 



'We shall replace —d^ by from now on and work in momentum space. 
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E.2 Calculation of the effective action 

We can now proceed in a very similar fashion as in section E.l. We can split the one loop effective 
potential in a few parts. A first part, ri^\ is the equivalent of (detP)"'/^ in expression (191) 



(1) . 



1 1 



3t)2ln 



M2 



/PP 



+(M2 + V^)2ln ^'+_V^ -2(M2 + 



PP' 



(203) 



The second part, the equivalent of (detX') ^Z^, is given by 



.(1) ^ 3(jV'-l) 
* 64ji2 



i-yi) 



--{m"- 2r ) - 7 (m^ - 2r ) +yf In 
e 6 /J 



+ 72 In 



i-yi) 



(204) 



whereby y i , and y3 are the solutions of the equation + {m^+ 2M'^)y^ + (X^ — pp"*" + 2M^np-)y+ 
M^X^ + l/2(p + p1")A,'* +M''m2 - m^ppt = and ^4 and ^5 of the equation / + 2M^y+M'^ - pp^ = 0. 
The third part is the constant term of the GZ action, 



From equation (138), we can calculate that^ 



(205) 



with 



1 + 



3 g'-N 1 
2 16ji2 e 



(206) 



so we find 
r(i) 



-d{N^ - 1)7^ = -4(iv2 - 1)/ - 4^ (AT^ - 1) g 1/ + ^ g/(^2 _ 1) 



X4 



^4 1 3 XV 2 X 



^ 'Ng^ ^ ^32Tfie 2 327t2 
The fourth part requires some calculation. We firstly find 



(207) 



4ZpZ2p \g^ g 



C3 CT4 



1 48(7V2- 1)2 

2 53iV 



53 1 TV _ pi \ PP 

6 el6ji2 ^ 24(Af2-l)2y g2 



(208) 



and secondly 



a' 



+ 



(J1O2 



4a'C'-3('2^2 ' 4ot/^/_^/2 g2 4a'^'_5('2 ^2 



2g- 



4ji2 



(209) 



so that 



^(i)_ 148(iV2-l) 



53iV 



53 1 Ng^ 



2 53 pi \ pp^ ^ow"* , OoM'' 



1 ° AT? 

6el67t2 * 24(Ar2-i)2y g2 2^2 g2 

/f4/Af2 _ 



I /UNCom'^ M4(Af2-l)2\ rim4 . , , 

+ - f + \,_7 )--^-M^ai+M2w^Xl- (210) 



967l2 



4;t2 



"For the explicit loop calculations of the Z-factors, we refer to [34]. 
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As a check on our results, we see that all the infinities cancel, so we find 



r(i): 



(Af2- 1 



--i6i,^[<*''-v'Pl''>''" 



2,„ M^-^ppt 



+ (M2 + Jppt)2in 



/ppT 



2 n2 



-2(iV2- 1) 



X4 3 A,' 

■ + 



4 



1 48(A^1) 

2 SSaT 



1-iV^ 



53 



Pi 



A'g2 2 323l2 
ppt 



(yv2 



24(Ar2-l)2; g2 



+ 1^-^ 7^ — - ^ - ATai +M2m2xi 



13 N 2g2 
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A' 



52 167t2 2 



(211) 
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